Modes of nonlinear acoustic transparency in the strained paramagnetic crystal 
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The propagation of the transverse-longitudinal acoustic pulses through a strained cubic crystal 
containing the resonant paramagnetic impurities with effective spin S = 1 is investigated. It is sup- 
posed that the pulses propagate under arbitrary angle with respect to the direction of the external 
static deformation parallel to the fourth-order symmetry axis. In this geometry, both the trans- 
verse and longitudinal components of the acoustic field have the high-frequency and zero- frequency 
spectral components. We show that a pulse can propagate in the modes different from the acoustic 
self-induced transparency. In particular, a pulse propagating in the mode of acoustic self-induced 
supertransparency substantially changes the populations of the spin sublevels, but its group veloc- 
ity remains almost equal to the linear velocity of the sound. If a pulse propagates in the acoustic 
extraordinary transparency mode, then its group velocity is substantially lower while the sublevel 
populations remain virtually invariant. Also, the modes of propagation under conditions of weakly 
excited spin transitions and large detuning of the pulse high-frequency components are identified. 

PACS numbers: 43.25.+y, 43.35.+d 



I. INTRODUCTION 

Historically, the coherent optical effects found their 
acoustic analogues after a while. Such a correspondence 
is easily traced, for instance, in the case of the self- 
induced transparency (SIT) phenomenon. Its discovery 
in 1967 and theoretical explanation in 1969 Q stimu- 
lated a search for the acoustic analogue. The experimen- 
tal observation at the liquid helium temperature of the 
acoustic self-induced transparency (ASIT) of the longitu- 
dinal hypersound in cubic crystal MgO containing para- 
magnetic ions Fe 2+ and Ni 2+ was carried out in 1970 by 
Shiren 0] , who gave also the theoretical consideration of 
this effect. Independently, ASIT was studied theoreti- 
cally for transverse hypersound in a system of paramag- 
netic impurities with effective spin S = 1/2 Soon 
after that, this effect was successfully revealed in a crys- 
tal LiNb0 3 alloyed by ions Fe 2+ @. 

Despite the existence of the optical-acoustic analogies 
mentioned above, the acoustic coherent effects do not ap- 
pear as the exact copies of optical ones. Obviously, the 
quantitative differences connected with a great distinc- 
tion of the light and sound velocities or the frequencies 
of the electric-dipole and spin-phonon transitions take 
place. It should be stressed that the qualitative differ- 
ences exist there. Thus, the optical wave is especially 
transverse whereas the acoustic one has the longitudinal- 
transverse structure. 

The interaction between transverse and longitudinal 
components of the acoustic pulse in a crystal without 
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anharmonicity of the lattice modes happens due to a 
presence of paramagnetic impurities. Its influence on the 
dynamics of the components is inessential, if the linear 
velocities of acoustic waves differ significantly, as it takes 
place in the most solids. However, the linear velocities 
of transverse and longitudinal sounds are close in the ion 
crystals of halogenide of alkaline metals Q- Since the 
components of the strain field interact effectively in this 
case, new features of their dynamics have to be expected. 
Indeed, new mode of the resonant acoustic transparency 
has been revealed theoretically in Ref. The propaga- 
tion of the transverse-longitudinal clastic pulses in this 
mode is accompanied by trapping of the population of 
the spin sublevels. Nevertheless, the deceleration in the 
velocity of the pulses is comparable to that in the case of 
ASIT of especially longitudinal or transverse hypersound 

si. 

The theoretical considerations of the nonlinear dynam- 
ics of transverse-longitudinal acoustic pulses were fulfilled 
basically for the paramagnetic ions with effective spin 
S = 1/2 IE 01 At the same time, it is well-known 
that the strongest interaction with the crystal lattice os- 
cillations takes place in the case S = 1 |f|. Since the 
spin-phonon interaction is greater here on two or three 
orders than for S = 1/2, the investigation of acoustic 
transparency of transverse-longitudinal pulses in a sys- 
tem of paramagnetic impurities with effective spin 5 = 1 
is more preferable from the point of view of experimental 
testing. 

The resonant modes of propagation of the transverse- 
longitudinal hypersound in direction parallel to external 
magnetic field and an axis of symmetry of the fourth or- 
der of the paramagnetic crystal with S = 1 have been 
studied in It was shown that these modes coincide 
with ones discussed in Ref. The distinctive feature of 
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the geometry considered there is that the roles of trans- 
verse and longitudinal components of the acoustic pulses 
are strictly various: the former component excites the 
spin-phonon transitions in the Zeeman triplet, the latter 
one causes the dynamic shifts of the frequencies of these 
transitions. 

Recently, it has been shown in the case of especially 
longitudinal acoustic pulses propagating under arbitrary 
angle with respect to an axis of the symmetry of the 
paramagnetic crystal, that both functions can be carried 
out by a single component of the strain field [ll|. I n 
this connection, the studying of the nonlinear modes of 
the acoustic transparency for the transverse-longitudinal 
pulses that propagate in a system of the effective spins 
5 = 1 under arbitrary direction to the symmetry axis is 
of great interest. The present report is devoted to allocat- 
ing these modes. Following to Ref. 0, we suppose here 
that the external magnetic field is absent, and splitting 
of the energy sublevels is created by the static deforma- 
tion of the paramagnetic crystal. For this reason, the 
degeneration of the spin states is removed incompletely 
in contrast to usual ASIT. 

The paper is organized as follows. In Section II, 
the system of material and wave equations describing 
the interaction in the strained crystal of the transverse- 
longitudinal elastic pulses with the paramagnetic impu- 
rities having effective spin 5 = 1 is derived. We use 
the semi-classical approach, i.e., the dynamics of the 
acoustic fields is assumed to obey the classical Hamil- 
ton equations for the continuous medium, but the spin 
subsystem is treated as a quantum object. The order 
of the derivatives in the wave equations on the strain 
field components is reduced with the help of the uni- 
directional propagation (UP) approximation. In Section 
III, the approximation of slowly varying envelopes (SVE) 
is applied to the particular case of the system obtained 
in the previous section. It is shown that the evolution 
of the acoustic pulses in the paramagnetic crystal is gov- 
erned by the well-studied equations of the long/short- 
wave coupling (LSWC) 0, ^1 • The pulse solutions de- 
creasing exponentially and rationally of the LSWC equa- 
tions are presented in Section IV. These solutions are ex- 
ploited to classify the modes of the acoustic transparency 
of transverse-longitudinal pulses in the deformed param- 
agnetic crystal. 



II. BASIC MODEL 

Consider a cubic crystal containing the paramagnetic 
impurities with effective spin 5=1. Let a transverse- 
longitudinal acoustic pulse propagate in the crystal along 
the z axis directed under arbitrary angle a with respect 
to one of its axes of symmetry of the fourth-order (z' 
axis; see Fig. 1). Also, we assume that an external field 
of static and uniform deformation employing the crystal 
is parallel to the z' axis. 

The Hamiltonian H of the spin-clastic interaction that 




FIG. 1: Scheme of the coordinate axes. 

arises in the most general case ^3 i s a function of the 
bilinear combination of spin operators Sj (j = x'^y, z'): 



(1) 



Coefficients Qij in this formula are supposed to depend 
on the components of the strain tensor 



i fdu m , du t 



2 \ dxi dx„ 

(m, I = x' , y, z'; U m are the components of the displace- 
ment vector U). In a representation of the eigenfunctions 
of operator 5 2 ', the spin matrices take the form [(| 
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The coefficients of the bilinear combination satisfy con- 
dition = Sij (8ij is the Kronecker symbol) if the ex- 
ternal deformation of the crystal is absent. Expanding 
these coefficients and function / in the power series on 
£ m i and retaining the first order terms, we obtain 



/ = f(S 2 ) + f'(S 2 ) 



(dQi 



\dS, 



n 1 1 



SiSjS m i , 



where subscript " 0" means that the values of derivatives 
correspond to nondeformed crystal. The first term de- 
pending on the Kazimir operator 5 2 = 5(5 + 1) = 2 is 
omitted in what follows as the constant addition to H. 
Then, Hamiltonian Q is expressed through the compo- 
nents of the spin-elastic interaction tensor 



G, 



ijlrn 



f'(S 2 ) 



dQi 



as given 



H — —Giji m £ m i(SiSj + SjSi) 
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In the case of the geometry we consider, the last equa- 
tion is rewritten in the following manner 

H = H Q + V, 



where Hq is the Hamiltonian of the effective spin in the 
7(0) 

'g'z' 



field £ 1,1, of the external deformation: 

1 



1 



Ho = G\\£ { °],Sl, = fej I 



(3) 



V is the Hamiltonian of interaction between the effective 
spin and strain field of the acoustic pulse propagating 
through the crystal 



V — G\\£ ZZ S Z 



Gx_ 
2 



£ xz {S x S z + S Z S X ) 



£yz(SyS z + S z Sy 



(4) 



Here we use notations G\\ = G 



z' z' z'z' 



G 



X'X'X'X' 



G 



yy , coo = G\\cy z , jh (H is the Planck constant) and 
G± = Gx'x'z'z 1 ~ G x ' x 'yy — Gyy Z > z >. Operators S x and 
S z are connected with S x i and S z i by means of the trans- 
formation of rotation 0, 



S x = S x ' cos a — S z > sin a, 
S z = S Z ' cos a + S X ' sin a. 



(5) 



For the consideration of the dynamics of the effective 
spins and the acoustic pulse to be self-consistent, we in- 
troduce the Hamiltonian of the elastic field 




(6) 



where p is the mean crystal density; p x , p y and p z are the 
components of momentum density of the local displace- 
ment of the crystal; an and a± are linear velocities of the 
longitudinal and transversal hypersounds, respectively. 
We assume in © that all dynamical displacements U x , 
U y and U z depend on variables z and t only. The integral 
is taken over the crystal volume. 

In accordance with the general scheme of the semi- 
classical approach 0, 0] , we describe evolution of the 
effective spin by the equation on the density matrix p : 



V,p\. 



(7) 



At the same time, dynamics of the acoustic field obey the 
classical Hamilton equations for the continuous medium: 



dp 

~dt = 
dU 
~dt 



S 

~sv 

6p~ 



- <V> 

<v>) , 



(8) 



<V>= n<V>dr, 



where 



n is the concentration of the paramagnetic impurities; 
< V> — Tr (pV) is the quantum average of V. 

Let the density matrix be represented in the next form 



P33 P32 P31 

P = I P23 P22 P21 

Pl3 PYX Pll 



(9) 



Eqs. J2J), and J3J) give us the following expression for 
the Hamiltonian of the spin-phonon interaction: 



(10) 
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where 



G\\ G 
Vn = —(1 + cos 2 a)£ zz ^- sin2a £ XZl 



V22 = G11 sin a £ zz + G± sin 2a £ xz , 



G11 sin 2a G± 
V12 = — y= — £ zz - -j={cos2a £ xz + icosa £ yz ), 



Gu 



sin 2 a £ zz + Gj_(sin2a £ xz /2 + i sina £ yz ). 



The physical mechanism of spin-elastic interaction in 
the case considered is the Van Vleck mechanism . The 
strain fields modulate the intracrystalline electric field 
in a location of the paramagnetic ions. Due to the 
quadrupole Stark effect, the static gradient of the elec- 
tric field causes a splitting of the quantum sublcvels of 
effective spin 5 = 1 that occur to be degenerate on the 
absolute value of its projection S z > (see © and Fig. 2). 
The components of the acoustic pulse excite the electro- 
quadrupole transitions between these sublevels and, as 
it follows from the expressions for Vn and V22, shift dy- 
namically the transition frequency. 

Using @) and © (|10fl . we deduce the wave equations 
on the strain field components: 



d 2 £ z 



nG\\ d 2 



dt 2 



11 dz 2 2p 8z 2 
sin 2a 



(3sin 2 a- 2)p 22 



+ 



V2 



d 2 £ 



12 2 
— a 



d 2 £ 



(P23 + P32 - Pl2 - P2l) 

+ sin 2 a(pi 3 -I- p 31 ) , (11) 

"Gl d 2 r . „ 

3sm2a p 22 



dt 2 x dz 2 4p dz 2 

+ \/2cos2a(p23 + P32 - P12 - P21 
+ sin 2a (p 13 + p 31 ) 



(12) 
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FIG. 2: Splitting of spin sublevels in the field of the static 
deformation (gray color), their position in a presence of the 
acoustic pulse and quantum transitions excited. The num- 
ber of the level and the corresponding value of projection S z i 
of effective spin on the axis of the external deformation are 
indicated from the left and right, respectively. 



end of Section IV that strict observance of this condition 
is not very important from physical point of view. The 
case of the acoustic pulse propagation through deformed 
paramagnetic crystal with effective spin S = 1, in which 
the linear velocities of the components differ essentially, 
has been studied in Ref. [Tj- 

As far as the linear velocities coincide, we are able to 
simplify Eqs. (|1 1|> — by reducing them in the order of 
derivative with the help of the UP approximation [l8j . 
Indeed, suppose that the concentration of the paramag- 
netic impurities in the crystal is small: 

_ n{Gl + Gl) 
and introduce new independent variables 



d 2 £ 



dt 2 



f) 2 F 



.nG ± d 2 



2p dz 2 



sin a (p 3 i - p i3 ) 



cos a 



{pi 



P32 — P21 



Equations J7J and (|1 1|> — form the self-consistent sys- 
tem describing the nonlinear dynamics of transverse- 
longitudinal acoustic pulses in the deformed crystal con- 
taining the paramagnetic impurities with effective spin 
5" = 1. Since its analysis is very complicated in the 
general case, we restrict subsequent consideration by a 
specific model. 

Let us assume that the linear velocities of both acoustic 
waves are equal: 

<2|| = aj_ = a. 

This restriction is rather artificial, since a|| > aj_ in 
solids. As noted in the previous section, it is fulfilled best 
of all in the ion crystals of halogcnidc of alkaline metals 
Q. One of the representatives of such the crystals is, 
for instance, NaBr Using typical experimental pa- 
rameters of medium and acoustic pulses, we show at the 



z 

T = t , ( = T1Z, 

a 

which arc usually referred to as the retarded time and 



P23) ■ (13) slow coordinate, respectively. Obviously, we have 



Of 



d_ 



d_ 

dz 



ld_ 

a dr 
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Then, neglecting the terms proportional to 7/ 2 , we ap- 
proximately write 



dz 2 



1 d 2 

a 2 dr 2 



2l _d 2 



a drd( 

in the left-hand side of Eqs. (fTT|) — (fT~Sfl and 



dz 2 



1 d 2 
a 2 dr 2 



in the right hand side. Integration on r of the wave equa- 
tions obtained by this way in the terms of new variables 
gives after returning to the initial ones: 



d£ zz 1 d£ z 



dz a dt 4V2hpa 3 



d£ xz , 1 d£ xz nG± 



i(pi2 - P21 + P23 - P32)(^o sin 2a - 2G±£ XZ ) + 2G± (V2 sin a {p 33 - P11) 

- COS a (p 12 + P21 + p23 + P32))£y 



dz a dt AV2hpa 3 
d£ yz 1 d£ yz nG± 



dz a dt 4:V2hpa 3 



i{pi2 - P21 + P23 - P32)(ft^o cos 2ci + G\\£ zz ) + Gj_(v2cosa (/9 33 - pn) 

+ Sina(pi2 + P21 + P23 + P32))£yz 

V2(pu - p 3 3)(G_L cosa£ KZ + G\\ sina£ zz ) + (p 12 + p 2 i + P23 + P32) 
x {Thuq cos a — G± sin a £ xz + Gu cosa£ zz ) 



(14) 



(15) 



(16) 
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It is assumed hereinafter that the components of the 
strain field and the density matrix elements satisfy fol- 
lowing conditions at t — ► — oo: 



J xz, yz, zz 



0, 



d£ 



xz, yz, zz 



0, Pjj 



Wi 



Pjk 







(j, k — 1, 2, 3; k ^ j). These conditions correspond to 
the pulsed mode of propagation of the sound in a crystal. 
Note that the evolution equation on density matrix 
is used to exclude the time derivatives of its elements in 
the right-hand side of Eqs. ifli jl -lfTB )! . 

The equations we came to are yet complicated for the 
rigorous analysis. Fortunately, they can be reduced in 
particular case to a system describing the interaction of 
the transverse-longitudinal acoustic pulses with two-level 
quantum particles. The investigation of this system will 
be carried out in the next sections. 



III. SYSTEM OF THE LSWC EQUATIONS 

Let the populations of the quantum states with S z > = 
±1 be equal: 

P33 = Pn- (17) 

This assumption is compatible with Eqs. (7J), (f 141) lllil) 

if following conditions are imposed 



P31 — Pl3, P31 



~Pl2, £y 



0. 



(18) 



It is checked immediately that conditions (|18|) are suffi- 
cient to keep Q17JI. Also, they give us relation 



P13 = Wi - /On. 



(19) 



Taking into account Eqs. I|17 p — 1|19| ) . we rewrite J7J and 
lO-JTU in the next form: 



dW 
~dt 



iEj_(a -a*), 



da 



dt 

dEj_ 1 dE ± 
— H 

dz a dt 



i(uj + 2E\\)a + 2iE x W, 
nA 



(20) 
(21) 



Ahpa 



■{lo + 2Eii+SE ± ) 



x(a-a*), (22) 



dEi 



II 



dz a dt 



1025,1 =i nA^ [s(u + 2 E\\+6E ± ) 



8hpa 3 



d 2 E ± ] {a -a*). (23) 



Here 



W= P22, cr = V2p21, 

Gj_ Gii 
E± = — cos 2a £ xz + — sin 2a £ zz , 
n 2n 



(24) 



G\ 

E\\ = -r cos 2a £ z 

2h 



G 



sin 2a £ xz 



and 



G\ cos^2a 



C' 2 
G ll 



sin 2a, 



S = 



2G\ 



^11 



2A 



sin 4a, 



V2 

d=^—G±_G 



II- 



System l|20ll - H23|) is remarkable from the point of view 
of its physical applications. In the case d = S = and 
En =0, it coincides with well-known reduced Maxwell- 
Bloch equations (RMB) for isotropic two-level medium 
[Isj . If d = 0, then we can put E\\ = -5Ej_/2, and the 
system appearing is nothing but the RMB equations for 
anisotropic medium [Tfij . which describe also the prop- 
agation of one-component acoustic pulses through the 
deformed paramagnetic crystals 0|- At last, Eqs. lf277jl ~ 
(|23|l are equivalent to the system of material and wave 
equations derived in Ref.0 under consideration of the dy- 
namics of transverse-longitudinal acoustic pulses in para- 
magnetic with effective spin S = 1/2 in the external mag- 
netic field presence. The complete investigation of this 
system turned out to be awkward. The simplest station- 
ary solution found corresponds to the case 5 — only. It 
describes the propagation of the extremely short acoustic 
pulse that possesses no well-defined carrier frequency in 
direction parallel to magnetic field. In contrast to this, 
the pulses containing of the higher-frequency components 
will be considered below. 

Owing to the last remark, further simplification of sys- 
tem (|2l))l (|2~3l will be achieved by applying the SVE ap- 
proximation. Comparing the right-hand side of Eqs. llL'L'li 
and l|23|) . we conclude that field En should have, in the 
general case, the high-frequency component proportional 
to Ej_. Taking this into account, we make use of the next 
representation 



a = R cxp[iuj(t — z/a)] , 
E± = (f2j_/2) exp[ioj(t — z/a)] + c.c, 
E ]] =(n ll -SE ± )/2, 



(25) 



where u> = luq — A is the carrier frequency (| A| -C ujq); R, 
and fin are slowly varying functions of t and z in the 
standard sense pfll ]. Substituting these expressions into 
Eqs. and disre sardine the nonresonant terms 

give us following equations 



dW 
~dt 



-(n* x R-tt x R*), 



— =i(A + flii)R + iQ x W, 
dt 



dn ± 1 dfl ± 



dz 

<9f2|| 
dz 



a dt 

1 0O|| 

a dt 



= -i(3±R, 
dW 



dt 



(26) 
(27) 
(28) 
(29) 
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where 



Pa 



najpA 

2hpa 3 ' ^ 



nAd 2 
Ahpa 3 



Note that we neglect the high-frequency addition to de- 
tuning A in the right-hand side of Ea. I|27[) . An influ- 
ence of this term on the dynamics of the pulses, whose 
duration is much greater than the oscillation period, is 
not significant because its average value over the pulse 
length tends to zero. In the case of the RMB equations 
for anisotropic media |19( . this fact was established in 
Ref.EB 

System (|26H - <|29[) we finally obtain differs only by the 
notations from the LSWC equations that have been in- 
vestigated in details due to their importance for theoret- 
ical study of the nonlinear dynamics of two-component 
electromagnetic pulses in the anisotropic resonant media 
[l2l IT3I ] . This reveals to us one more correspondence be- 
tween coherent optical and acoustic phenomena, which 
will be used in the next section for classification of the 
modes of the acoustic pulse propagation through strained 
paramagnetic crystal. It should be noted here that the 
system of LSWC equations, as well as its gauge equiva- 
lents, arises in various physical problems (see, e.g., refer- 
ences in Ref.lH. These numerous applications point to 
the universal physical character of the LSWC system. 

One of the features of Eqs. lj2T)j) - lj2l?|) is that the field 
variables play in them totally different roles. Namely, 
Q_l causes the quantum transitions, whereas Qj_ shifts 
dynamically the transition frequency. It is easy to see 
that there exists relation between the field components: 



fi N = - — \n ± \ 2 +F{t-z/a). 



4wi 



(30) 



The last term on the right-hand side generates the phase 
modulation of R and il± only and vanishes after appro- 
priate change of the variables. Without loss of a gener- 
ality, we put F(t — z/a) = in what follows. 

At the end of this section, we find the expressions for 
transverse and longitudinal components of the strain field 
through the variables of the LSWC equations. It follows 
from Eqs. ijHJ and that 



reveal us an asymmetry on the polarity of acoustic sig- 
nal: the signs of the zero harmonics of its transverse and 
longitudinal components arc determined by a and the 
type of the external action (tension or compression) on 
the crystal. It is remarkable that the ratio of the ampli- 
tudes of the zero harmonic of £ zz and £ xz : 



G ± 



G 



cot 2a, 



II 



is equal to inverse ratio of the first harmonic case. This 
fact can be used in measuring the constants of the spin- 
phonon interaction in the paramagnetic crystals. 



IV. THE MODES OF ACOUSTIC 
TRANSPARENCY 

The solution of Eqs. Ij26ll Il29|) that describes the prop- 
agation of the transverse-longitudinal acoustic pulse is 
written as given (see Ref. 0) : 



Here 



M 



± = VM exp(i$), 



d 2 

flu = M, 

4w 



89 



t 2 [g - a + v /l + (.g-a) 2 cosh2C 



Wqz — arctan 



tanh ( 



const, 



(33) 



(34) 



(35) 



(36) 



(37) 



a = At p , W 



1-3W 2 



f 



f — 



cos2a^ilj^ cxp[iuj(t — z/a)] + c.c.^j 



2G± 



sin2a On, (31) 



sin 2a(p,± exp[iw(i — z/a)] + c.c^j 



— — cos 2a Q 



G 



II 



II- 



(32) 



2bJ Tp 



g - a + V 1 + (3- a) 2 , 



v„ =o(l 



a(3±r 2 
1 + a 2 



Wn 



(38) 



These formulas display that both the components of the 
acoustic pulse have high-frequency and zero-frequency 
components in the general case. Also, Eqs. (I3UII 1321) 



The free parameters of the pulse presented are A and t p . 
For the sake of convenience wc suppose hereinafter that 

T p > 0. 
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As it follows from Eqs. IM-il) and (|37|) . component fix 
has the phase modulation, which leads to a local non- 
linear chirping of the carrier frequency, to — > u>i oc = 
lu + StOnon, where 



UUJ non — 

OT 4 



(39) 



Taking into account the dynamic shift of the transition 



frequency luq 



1 J 



0|| (see Eq. J2Z|), we come 



to the following expression for effective detuning A e j 
component fl± from the resonance: 



of 



A- 



I n ll- 



(40) 



Defining the pulse length T p as the double deviation 
from the zero point of t — z/v g , at which jfijj is half its 
maximum value, we obtain from formula (|33() : 



T p = t p arccosh 4 + 3 



g- a 



(41) 



In the SVE approximation, both pulse length and nonlin- 
ear shift of the carrier frequency must obviously satisfy 
the conditions ojqT p 1 and |£u; non | <C u>o- It can easily 
be shown that these inequalities are valid if ujqt p 3> 1 
and loqt p a — g. The last condition is necessary in the 
case a — g 3> 1, evidently, and can be fulfilled only when 
d 2 » 1. 

Considering the limit t p — > oo in formulas (|33|I - H38(I 
gives us rationally decreasing pulse solution: 



Let us begin with the modes, in which the pulse excites 
the paramagnetic impurities strongly. Since W 2 + \R\ 2 
is independent of t, the strongest degree of excitation 
happens in the case, when the value of variable W in the 
pulse center differs by the sign from its value at absence of 
the pulse. By virtue of Eqs. (|35|l and (|36J) . this condition 
can be formulated as 



r p 2 |lMC = 0)| 2 = 4(l + a 2 ). 



Performing some simple algebra, we put it into the next 
form 



(a 2 - 2ga+l) 2 = 0, 



which yields 



2 \ a 



(44) 



This condition can take place if a > only. It can also 
be rewritten as a relation between the pulse detuning A 
and parameter t p : 



4^o 
d 2 



A - A 2 



(45) 



d 2 {l + in 2 £) 



exp^xWbz/A), 



W=[l 



8k 2 



(i + K 2 ) 2 (i + K*e) 



(42) 



(43) 



where 




Parameter k in Eqs. I|42|l and (|43|l is supposed to be real. 
This imposes a constraint on detuning A of the rationally 
decreasing pulses: A > A r , where A r = 2ujo/d 2 . 

Now, we use the expressions presented above to 
identify the modes of the propagation of transverse- 
longitudinal acoustic pulses in deformed paramagnetic 
crystal. These modes resemble the ones studied in Ref.fl3l 
for the case of two-component electromagnetic pulses. In 
what follows, it is assumed for the sake of concreteness 
that the paramagnetic impurities arc in thermodynamic 
equilibrium prior to the pulse passage (—1/2 < Wo < 0). 



Setting the right-hand side to zero, we find the interval 
of admissible values of the detuning: < A < A rn = 
2A r . The curve of the strong excitation (|45|) is plotted 
in Fig. 3. 













V 1 








in 






IV 






V 




i 







2A r 



FIG. 3: Curve of the strong excitation. Domains of the 
pulse parameters corresponding to different modes of acous- 
tic transparency: (I) ASIT; (II) ASIST; (III) AEOT; (IV) 
ANNT; (V) APNT. 



It is seen from Eq. Ij43(l that the strongest excitation 
in the case of rationally decreasing pulse takes place if 
K = 1. Obviously, this agrees with l|45|) . 
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Substitution of relation $Q into Eqs. l|3"3" )l -l|3"8" |) yields 
|fi x | = - 



2Vl + Q' 2 



f2n 



W=[l 



T pV /l + (l + Q' 2 )sinh 2 C 
4A 

1 + (1 + a 2 )sinh 2 C' 
2 



l + (l + a 2 )sinh 2 C 



(46) 



Wo, 



The corresponding pulse length is 



T p = 2t p arcsinh 



(47) 



Accordingly to Eqs. (|40|l and (|46|l . the effective detuning 
of fij_ from the resonance equals A at the edges of the 
pulse and — 2A at its center. Thus, we can say that 
this component is resonant with quantum transitions on 
average over the pulse length. As a result, the largest 
possible change in the population of the spin sublcvels is 
achieved. 

If the detuning is very small (a<l,j>l), then Eqs. 
11401 can be represented as 



|0_l| = — sechC, fin = -4Asech 2 C 



II 

W=(l - 2sech 2 C) W , 
-l 



(48) 



o(l - a(3 ± r^W 



One can see that |Oii| <C |f2j_|, and the phase-modulation 
depth for component fl± is much smaller than its input 
spectral width: |<5cj, io „| <C 1/t p . In the exact resonance 
case (a = 0, g — > oo), these relations give us well-known 
expressions of the ASIT theory 

Formulas (|46f) and (|39[) show that an increase of de- 
tuning (and therefore a) leads to larger amplitudes of 
the pulse components and to a deeper chirping of f2j_ to- 
ward lower frequencies. Since the pulse length decreases, 
the profiles of both components become sharper. 

Also, it is necessary to note that the group velocity 
approaches linear velocity a of acoustic wave as detun- 
ing increases. Nevertheless, the strong excitation of the 
paramagnetic impurities takes place: the largest possible 
change in the level population is reached at the center of 
the pulse. This is evidently explained by the fact that a 
growth in the amplitude of pulse increases its power. The 
ensuing higher rate of excitation/de-excitation processes 
leads to a higher soliton propagation velocity. 

It is clear from l|44|) that each particular value of g 
corresponds to two distinct values of a. We associate 
the domain of relatively small detuning (a < 1) with the 
ASIT mode, since it is implemented in the case a = 0. 
When a > 1, we say that the pulse solutions given by 



(|46|l propagate in the mode of acoustic self-induced su- 
pertransparency (ASIST), thus emphasizing the fact that 
the group velocity decrease is lower as compared to the 
ASIT mode while excitation is equally strong. Since 
|A e /| -C in the SVE approximation, the ASIST mode 
must be most strongly manifested if d 2 3> 1. The do- 
mains of the values of the pulse parameters, which cor- 
respond to these modes, are schematized in Fig. 3. 

When j>1 and a ^> 1, the amplitudes of components 
fl± and fin of the ASIST pulse and the phase-modulation 
depth for Qj_ reach their limits 2A m , 4A m and A m , re- 
spectively. In this case, (|47H gives us the following esti- 
mate for the pulse length: 



oj ' 



This time scale corresponds to the time scale of phase- 
modulation localization. A deeper phase modulation 
combined with a decrease in the corresponding localiza- 
tion time scale and a shorter pulse length can be inter- 
preted as an effect of the spectral supcrcontinuum gen- 
eration. Indeed, the Fourier transform of Q33|l and l|34(l 
defined as 



F±,\\(y) 



dt 



yields 



\F ± (y)\=2v 



V9 



1t p v/2 



\/l + (g-a) 2 cosh(7rr p ^/2) ' 



\F\\{v)\=to 



sinh(0T p ^/2) 
sinh(7TTpi//2) ' 



where 9 = arccot((7 — a) (0 < 8 < ir). The absolute 
values of the Fourier transforms of the components fi|| 
and Qj_ reach their maximum values at v — and v = uq, 
respectively, where 



2 , 7T + I 

In 

TTTp 7T — I 



If a — g 3> 1 or, equivalently, — > n, then the spectral 
width of the pulse is Suj ~ — 0)t p and the max- 

imum of the spectral energy distribution of component 
Q± is reached at frequency lj — uq, which is much lower 
than the carrier frequency lj. In this case, despite a large 
linear detuning from resonance (r p A 3> 1), a substantial 
nonlinear spectral broadening (t p 5lu 3> 1) and dynamical 
shift in the transition frequency luq — > ufj lead to a gen- 
eration of the resonant Fourier components (photons), 
which stimulate the quantum transitions. 

Let us now consider the modes of the acoustic pulse 
propagation in the case g <^ 1. Although the paramag- 
netic impurities arc excited weakly under this condition, 
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there exist remarkable features in their interaction with 
the acoustic pulses. Since lu t p ^> 1 in the SVE ap- 
proximation, the value of \d\ must be sufficiently large to 
ensure that ujQT p /d 2 <C 1. This can be fulfilled for some 
values of a if |Gj_| 3> |G||| or vice versa. 

First, we assume that the detuning from resonance is 
large (\a\ > 1). If a < (or u > lj ), then Eqs. 
P)l give 



n 




ii 



sechC, 



sech C, 



(49) 



2<? 



W = [ 1 - sech 2 C ) Wo, 



Comparing these expressions with (|48|l . we see that the 
amplitude of Qj_ is much smaller than in the ASIT mode, 
whereas the respective f2||'s are comparable. We also 
note that 

l«x(C = 0)| rr-T 



|n N (C = o)| vai " 

i.e., the amplitude ratio can have an arbitrary value. The 
paramagnetic impurities remain almost unexcited as the 
soliton described by 1(49(1 propagates through the crystal, 
and the soliton velocity decreases only very slightly. The 
phase-modulation depth is also small (|<5u> ncm | <C 1/t p ), 
and the effective detuning A e / only increases as the com- 
ponent f2|| is generated. This leads to an even weaker 
excitation of the impurities as compared to that induced 
by the input pulse. 

Now, let us suppose that j « 1 and a ^> 1. In this 
case Eqs. p5 )l -(|5g )l imply 



IfiJ = 



n,i = - 



t v \J\ + 4a 2 sinh 2 C' 
8A 



W = 1 



1 + 4a 2 sinh 2 C : 

8g 

a(l + 4a 2 sinh 2 Q 
w„ = v r . 



(50) 



Wo 



Note that expressions ((46(1 and (| 5 Of) are somewhat sim- 
ilar: in both cases, the solitons are sharply peaked, and 
their propagation velocities are nearly equal to a. How- 
ever, these modes are essentially different in terms of be- 
havior of the paramagnetic impurities. Whereas they are 
strongly excited as the soliton described by ((46(1 propa- 
gates through crystal, no significant excitation is caused 
in the case of l(50"}) . Indeed, since the effective detuning 
A e f of the pulse described by Eqs. (|50"|l is — 5 A at its 
center, it is not resonant with the paramagnetic impuri- 
ties on average over the pulse length. However, according 



to (14911 and 1150(1 . the excitation of the paramagnetic im- 
purities at u> < ujq being relatively weak, is still stronger 
than that at lj > u>q. The reason is that the effective de- 
tuning decreases toward the pulse center when lu < u)q, 
owing to the component On , and increases when u> > luq. 
Thus, a comparison of ((49(1 with 1(50(1 demonstrates ob- 
vious asymmetry with respect to detuning of f2j_ ■ Since 
a < for the pulses described by l(49|) and a > for ones 
described by (|5TJ|) , we refer to the corresponding modes as 
acoustic negative and positive nonresonant transparency 
(ANNT and APNT), respectively. 

If g <C 1 and \a\ <C 1 (detuning is small), then Eqs. 
((33(1 138|) lead to expressions identical to those found in 
Ref.[l2: 



|fix| = -- 



/ 2 ff sech 1/2 2C 



flu = sech2£, 

W = (1 -4 5 sech2C)Wo. 



(51) 



The expression for group velocity coincides with one 
given in Eqs. (|4*5|) corresponding to the ASIT mode. 

Here, we have |0||/J7_l| 2 ^> 1 in the center of the pulse. 
The paramagnetic impurities are not excited, since the 
effective detuning is large, 3/r p . However, the propa- 
gation velocity decreases as in the case of strong exci- 
tation at a ~ 1. Such a deceleration of the pulse is 
explained by the dispersion properties of medium within 
the higher-frequency component bandwidth. According 
to ^3 ) where this effect was studied in details for the op- 
tical solitons, we call the mode considered an acoustic ex- 
traordinary transparency (AEOT). The roles of the ordi- 
nary and extraordinary pulse components are played here 
by Qj_ and f2||, respectively. The existence of this mode 
for transverse-longitudinal acoustic pulses has been re- 
vealed also in the case, when the splitting of the spin sub- 
levels is produced by the external magnetic field @, El • 
The domains of existence of the modes with g <§C 1 are 
presented in Fig. 3. 

It should be noted that expressions (13311 138|) do not 
change their form in a presence of detuning between <Z|| 
and aj_. The influence of this detuning is inessential if 
condition 



1 

aj_ 



1 

"II 



< 1 



is valid. Substituting the definitions for (3± and A into 
(138(1 . we obtain 



a/3j_T 2 \a±_ 



1 

"II 



hpa?. 



nGf,uj T 2 \aj_ 



Taking for Fc 2+ :MgO 0, IE O 113, n ~ 10 



17 , 



L0 



10 



10 , 



10 



-13 



erg, p 



2g/cm , a 



II 



5-10 5 cm/s, ffl||/ax ~ 1-5 and r p ~ 10~ 8 s, we find e ~ 0.1. 
This estimation shows that we can neglect detuning be- 
tween the linear velocities of transverse and longitudinal 
acoustic waves for typical values of the ASIT parameters. 
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TABLE I: Characteristics of the modes of acoustic trans- 
parency 



Mode 




\a\ 


9 


Vg 


|OlA2||| 


Excita- 
tion 


ASIT 


-2A 


« 1 


> 1 


VASIT 


> 1 


strong 


ASIST 


-2A 


> 1 


> 1 


< a 


~ 1 


strong 


ANNT 


w A 


> 1 


< 1 


a 


arbitrary 


weak 


APNT 


-5A 


> 1 


< 1 


a 


< 1 


weak 


AEOT 


»|A| 


< 1 


« 1 


> VASIT 


« 1 


weak 



"Note: A e f is the effective detuning of the high-frequency com- 
ponents at the center of acoustic pulse; A e j is compared with A 
in the AEOT mode on absolute value. 



V. CONCLUSION 

In this paper, we have investigated the soliton modes 
of the acoustic transparency in a strained cubic crystal 
containing the resonant paramagnetic impurities with ef- 
fective spin S — 1. It is supposed that the linear veloc- 
ities of transverse and longitudinal sound are close, and 
the pulses propagate through the crystal under arbitrary 
angle with respect to the direction of external deforma- 
tion parallel to the fourth-order symmetry axis. 

We have allocated five modes of acoustic transparency, 
which differ by the propagation velocity of the transverse- 
longitudinal pulses and degree of excitation of the param- 
agnetic impurities (see Table I) . The acoustic self- induced 
transparency is characterized by strong excitation and 
substantial deceleration in the pulse propagation veloc- 
ity relative to linear velocities. Self- induced supertrans- 
parency differs from ASIT in that the decrease of veloc- 
ity is small, but the paramagnetic impurities are strongly 
excited as well. The solitons of ASIST have larger am- 
plitudes and smaller lengths as compared to ASIT soli- 
tons, and their high-frequency spectral components are 
strongly modulated. The carrier frequency of the trans- 



verse and longitudinal components in this mode is lower 
than the resonant frequency. The modes, in which the 
trapping of the populations of the spin sublevels takes 
place, are also identified. The pulses propagating in the 
acoustic extraordinary transparency mode are character- 
ized by small detuning of the high-frequency components 
and dominant role of the zero-frequency ones. Their 
group velocity substantially changes and may become 
comparable to that of pulses in ASIT and ASIST modes. 
In the acoustic positive and negative nonresonant trans- 
parency modes, the pulse velocity changes insignificantly, 
and the absolute value of detuning is large. The most 
substantial difference between these modes concerns the 
behavior of the effective detuning of the high-frequency 
components. In the ANNT mode, it remains virtually 
constant. If a pulse propagates in the APNT mode, then 
the effective detuning changes sign due to an influence of 
the zero- frequency component. Since the local frequency 
passes through a resonance, a slightly stronger excitation 
occurs in this case, and the pulses are sharply peaked, as 
in the ASIST mode. 

In this study, we ignored inhomogencous broadening 
of the spin sublevels. An investigation allowing for this 
effect may throw light on the pulse area theorem, provid- 
ing a basis for analysis of the stability of the transverse- 
longitudinal acoustic pulses. Also, it would be interest- 
ing to identify the distinctive features of the dynamics 
of picosecond transverse-longitudinal acoustic pulses in a 
system of effective spins 5 = 1. Unlike the present con- 
sideration, the SVE approximation is inapplicable to this 
case. 
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